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A multiple scattering formulation for the electromigration wind force on atoms in dilute 
alloys is developed. The theory describes electromigration via a vacancy mechanism. The 
method is used to calculate the wind valence for electromigration in various host metals 
having a close-packed lattice structure, namely aluminum, the noble metals copper, silver 
and gold and the Ad transition metals. The self-electromigration results for aluminum and 
the noble metals compare well with experimental data. For the 4c? metals small wind valences 
are found, which make these metals attractive candidates for the experimental study of the 
direct valence. 



I. INTRODUCTION 



An electric field, applied to a metal sample, causes an atomic current, besides the common electric current. This 
process is called electromigration and has been studied for many years. Already in 1930 Coehn and Spechtu reported 
a flow of hydrogen atoms towards the cathode in palladium. This migration of interstitial atoms is a rather harmless 
manifestation of the effect. It offers the possibility to influence the concentration profile of impurity atoms along a 
sample. However, a net flow of host atoms can also be the consequence of applying an electric field. Clearly, this 
reduces the conductivity of a metal wire because of the formation of voids on one side and hillocks on the other side. 
In wires of normal size this turns out not to be of great importance, because the atoms do not move very fast. In thin 
films, however, the situation is different. Thin films can carry higher current densities, because the heat, produced 
by the current, is transfered more easily to the environment. Therefore they can endure a current density greater 
than 10 7 A/ cm 2 , while a bulk sample would melt when the current density is 10 4 A/cm 2 . Electromigration damage 
in aluminum films is a well-known example. 

The migration of atoms turns out to be influenced by the presence of impurities. Most important is the effect of 
copper doping of Al films, which greatly reduces the electromigration damage. This effect is known for a long time 
without understanding it, so impurities have been added using the method of trial and error. Furthermore, although 
the damage appears too fast for a device user, it is still not fast enough for an experimental physicist. So, the effect 
is studied under so-called accelerated conditions, high temperatures and high current densities. Unfortunately, the 
extrapolation to user conditions is not trivial, because the total atomic flow is a result of a number of competing con- 
tributions with different activation energies. All of this complicates the picture and cries for a theoretical framework, 
to which this article is aimed to contribute. 

Two ingredients are essential in the process of electromigration, diffusion and a driving force. The motion of 
atoms due to diffusion is random. At lower temperatures diffusion along grain boundaries is dominant, but at high 
temperatures bulk diffusion becomes the most important contribution. For accelerated conditions it is the leading 
mechanism and the atoms move by exchanging positions with vacancies. The second ingredient, the driving force, 
introduces a bias into the random motion of atoms. One contribution to that force is calculated in the present work, 
namely the wind force due to scattering of electrons by the moving atom, which can be either_a host atom or an 
impurity atom. A well-established quantum-mechanical expression for the wind force is available!! Besides the wind 
force, the driving force has a direct part. If the nuclear charge is not completely screened, the electric field-directly 
pushes the atom towards the cathode. This force has been the subject of a controversy for many yearsBa. Both 
contributions are proportional to the electric field and the total driving force can therefore be written as 

F = F win d + Fdiroct = (Zwind + Zdiroct)eE = Z*e~E. (1) 

The effective valence Z*, which is the sum of the wind and direct valence, has been measured for a lot of systemsi. 
Often the wind force dominates the direct force and, depending on the system, the resulting total force can push the 
atom either to the anode or to the cathode. 

The calculation of the wind force requires knowledge of the electronic structure of the alloy. An impurity in a 
dilute alloy, in which the concentration of impurities is low, only interacts with its direct environment of host atoms. 
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The electrons feel a potential, which only differs from the host potential in a small cluster of atoms containing the 
migrating atom, the vacancy and the surrounding host atoms affected by charge transfer and lattice deformation. A 
Green's function formulation of multiple scattering theory is used to calculate the electron wave function in the dilute 
alloy. This formalism-is an extension of the Korringa-Kohn-Rostoker (KKR) method for the calculation of the band 
structure of a metalBQ. The alloy is described with respect to a reference system. The conventional choice of this 
reference system is the host metal, but then the electronic structure of the important saddle point configuration in 
electromigration cannot be described. This is due to the fact, that the atoms in the alloy system cannot be mapped 
onto the ones in the host system one-to-one. As proposed by LodderQ an intermediate void system is used, which does 
not contain any scatterers in the perturbed region. The original formulationO, in which the cluster of perturbed atoms 
is considered as a whole, was recently applied in a preliminary study of electromigration in copper and aluminumO. 
This formulation is rather cumbersome and in the mean time we found a simpler single site formulation, in which 
each atom in the cluster-is considered .separately, analogous to the formalism used for the description of interstitial 
impurities by OppencerEij and van Eklill. The validity of this method is not limited to the particular problem of 
substitutional electromigration. In all problems concerning dilute alloys the wave function can be calculated this way. 
The multiple scattering expression for the wind force is elaborated in section [n]. This section contains two sub- 



sections, devoted to the wave function and the Green's function matrix respectively. In section III some details of 



the calculation are given. The formalism is applied to metals with a close-packe d latt ice structure in section IV 



Self-electromigration and impurity electromigration in Al are t reate d in subsection IV A , self-electromigration in the 



noble metals as well as impurity migration in Ag in subsection [V B and self-electromigration in 4<i transition metals 
in subsection IV C. In section ^ the main results of this paper are summarized. 



Throughout the article atomic units are used, such that Ti = 2m = 1. Exceptions are stated explicitly. 

II. MULTIPLE-SCATTERING EXPRESSION FOR THE WIND FORCE 

In this section the quantum-mechanical equation for the wind force of Sorbello et al.! 

F win d = < **l - Vr,^** >. ( 2 ) 

ft 

which traces back to the pioneering work of Bosvieux and FriedelB, and in which 

5f(k) =er fc E-v fe d/ (e fc )/de fe , (3) 

is expressed in computable quantities. In Eq. (^|) the alloy electron wave function, is supposed to be constructed 
from the corresponding Bloch function, labeled by crystal momentum k and band index n, combined in k = (k, n), 
and having an energy eigenvalue eft. In a simple ballistic picture momentum is transfered from electrons to the 
atom because of scattering, resulting in the wind force. The momentum transfered to this atom, with a potential v p 
centered at position R p in the alloy, is proportional to the expectation value of the force operator -V^tip. When 
no electric field is present, the electrons are distributed according to the Fermi-Dirac function /o(efc), leading to an 
exact cancellation of forces. In other words, each push from the electrons in one direction is compensated, on the 
average, by one in the opposite direction. This symmetry is broken by the presence of the electric field, which alters 
the distribution function by Sf(k). The first and leading term of this deviation, Eq. (^), is linear in the applied 
electric field E and proportional to the electron transport relaxation time r^. In Eq. (^|) e is the elementary charge 
and Vfc is the velocity of the electron. As usual, in all applications, the transport relaxation time is taken constant 
throughout the Brillouin zone. All other quantities can be extracted fronipthe electronic structure. The method for 
the calculation of the electron wave function originally provided by LodderEl will be used in a modified way. 

In the present formulation the crystal is divided into cells. For the host crystal these cells are the well-known 
Wigner-Seitz cells containing local potentials w h (x). For the alloy these cells may have a different shape in the region 
of the impurity cluster and they contain local potentials v p (x) . The local potentials v h and v p consist of a spherical 
atomic potential surrounded by an interstitial region, where the potential is constant. In order to evaluate the matrix 
element appearing in Eq. (^), the wave function is expanded in the p-th cell, centered at R p . In such a cell the 
wave function 'J'fe(r) = ^ / / £ (x + R p ) can be written as a linear combination of regular solutions R P L of the Schrodinger 
equation with local potential v p , so 

vK fe (x + R p ) = ^c fcpi i?P(x), (4) 
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where L = (I, m) combines th e an gular momentum and magnetic quantum numbers. The expansion coefficients Ck P L 
will be derived in Subsection II A . The basis functions R P L are constructed such that outside the atomic sphere they 
have the free space form 

i^(x)=i £ (x)-i^;tj i ,fe+(x), (5) 

L' 

where Jl(x) is the product of a spherical Bessel function ji(nx) and a spherical harmonic Yl(x), ftj(x) is a similar 
product for the spherical Hankel function hf(nx), k — \f~E and t p is the scattering matrix for the potential v p embed- 
ded in free space. For a spherically symmetrical scatterer this matrix is diagonal and equals t p LL , = — sin rjf e lT,v i Sll' , 
the r) being the phase-shifts. 

In that case the wind force expression (|J) can be elaborated to the form 

l+l 

F wind = ^/(fc)2Re^ J2 4 P LD L;i+1 , mi sin(7 7i+1 -7 7i )e 4( " !+1 - r ' i) Cfcpi+ i, mi , (6) 

fe L mi=— (l+l) 

according to Eqs. (26) and (31) in Ref.i, in which the vectorial matrix D is defined as 

Dll' = / dxY L {x)xY L ,{x). (7) 



A. The wave function coefficients 

The alloy wave function coefficients that appear in Eq. (|]) will be expressed in terms of wave function coefficients 
c fcjL i which appear in the local expansion of Bloch function vf£ ost (r) = ^)j ost (x + Rj)> 

^ ost (x + R J -) = ^4°lX(x) ! (8) 

L 

where Ri are the regular solutions of the Schrodinger equation within the Wigner-Seitz cell labeled by j and contain 
t h instead of t p in their asymptotic form (see Eq. (^)). Unperturbed atomic positions Rj in the host do not necessarily 
coincide with one of the atomic positions R p in the alloy. The number of alloy sites may even differ from the number 
of host lattice sites. Therefore, although the alloy is described with respect to the host, the host system cannot be 
used as a reference system straightforwardly. An intermediate system is required, which will be referred to as the void 
system. This system consists of a void region, where the potential is constant, surrounded by host atom potentials. 
The void region is chosen to be as extended as the perturbed region in the dilute alloy. This system can serve as a 
reference system for both the host and the alloy, as it is represented symbolically in Fig. [I]. The relation between the 
host and alloy wave function is given via the void wave,function. The alloy wave function is expressed in terms of the 
void wave function by the Green's function expressionla, 

* fe = *r id + G void vP ^ k - ( 9 ) 

p 

Analogously to Eqs. (||) and (H) the void wave function can be written in terms of local basis functions, being spherical 
Bessel functions because of the constant potential in the void region 

*r d (x+R P ) = ^ c ^ L d j L (x). (io) 

L 

Inside the void, the void Green's function can be written as 

G void (x + R ?) ,x' + R P = - JK ^j L (x < )^+(x > )V 

L 

+5>(x)sr^ w 'i i (xo. (ii) 

LL' 
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This expression with the void Green's function matrix g vmd <PP is a straightforward generalization of the expression 
for the Green's function of free space 



G°(x + R P! x' + R P = -m^i i (x < )/i+(x > )5 p 

L 



LL' 



in which the free space propagation matrix B appears, defined as 

= B LL ,(R PP ,) = W- ; '-V(l - 5 pp/ )J2^"CLL'L"h+,(R pp ,). 



(12) 



(13) 



The void Green's function matrix will be given in Subsection [I B . 

Combining Eqs. (^|), (10) and (|TT| ) the following matrix equation can be found 

void \ "V-i /->void . \pp 
C kpL = Z^ 1 ~ 9 tfLL'Ckp'L'- 



(14) 



Details of the derivation will be given in Appendix [A|. Following the same procedure, void wave function coefficients 
carrying the host position label j can be expressed in terms of host coefficients cjyj 



c ^ = E( 1 -^° idt ' l )ii' c ^'- 

j'L> 



Since c v h °^ and are related by 



kjL 



void _ \ jpj void 
Li 



(15) 



(16) 



Eqs. (|14| ) and (15) can be combined to 



p'L' 



3'L'Li 



(17) 



It is noteworthy that in this equation two types of summations over the angular momentum occur. The first type has 
a natural cut-off because of the multiplication by i-matrices, and the second type, such as the summation over L\, 
in principle runs to infinity. The latter summation can be treated analytically, to be shown in Appendix |^, and the 
final equation used in the actual calculations, 



EC 1 - G void t) pp 'c kp , = 



p' 



6(k,R w )e 4k R «i h c)!° st - ^£ void < ro 't h c£°, st , 



(18) 



j' 



only contains summations with a natural cut-off. Angular momentum labels have been dropped in order to simplify 
the notation. The matrix 6(k, R) is defined by 



6(k,R) = Y,B(R-B. j ,)e- ik < R - n i'\ 



(19) 



i' 

in which the summation runs over all lattice sites. 

The host wave function coefficients occurring in Eq. ([Lq) are calculated by the KKR method according to 



£M 1£ ,(k)#c&£ = 0. 



(20) 
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This equation contains the KKR matrix 



M(k)=t fc -6(k), (21) 

in which 6(k) = 6(k, 0), defined by Eq. ©. „ 

In its original form the KKR method applies to muffin-tin potentials. However, as shown by NesbetEJ, it is possible 
to circumvent this restriction on the crystal potential in the same way as has been done in the derivation of Eq. . 
Now we will use this procedure in the calculation of the void Green's function. 

B. The void Green's function 

The void Green's function and the host Green's function G h are related by the Lippmann-Schwinger equation 

G void _ G h _ G h v *G void . (22) 

j 

The summation runs over the positions in the void. This equation is again evaluated in the position representation. 
The host Green's function is given by 

G h (x + R, , x' + R y ) = ]T i?L (x< )S L (x>)^ 

L 

+ Y,RlW l {,R l (*'). (23) 

LL' 

Outside the atomic sphere the regular solutions Rl of the local Schrodinger equation have the asymptotic form of 
Eq. (||) with t p replaced by t h and the singular solutions Sl equal —iKh\. First the void Green's function matrix 
£ void j j ^ carr yi n g host position labels j and j', will be expressed in terms of the host Green's function matrix QH , 
to be determined later. To that end Eq. ( p2| ) is elaborated at the positions x + Rj and x' + Rj/ , where both x and 
x' are located close to their cell centers Rj and Ry respectively. Following similar steps as in Appendix ^ one finds 
the relation 

gvoid, jj' _ gjj' _ gjjt ^.hgvoidjif ^24) 
h 

where the angular momentum labels have been dropped. The matrix t h is the same at all lattice sites, so the site 
label has been dropped. Iteration of this equation leads to 

gvoid ji' = gjj' _ gm (j*- 1 + g) h ) 2 g^'. (25) 

The host Green's function matrix g is calculated using the Lippmann-Schwinger equation 

G h = G° + G°^w 3 h G h , (26) 

j 

G° being the free-space Green's function. The summation runs over all host lattice sites. Working in the position 
representation and using Eqs. ( p"2| ) and ( p3| ) one finds the following matrix equation 

gjj' = B jj' + B jji t bgjij' _ (27) 
Taking the Fourier transform of this equation one finds 

g{k) = J2 ^'"''' kH ' = 6(k) + b(k)t h g(k). (28) 

j' 

This equation can be solved for g (k) 



5 



G(k) = (1 - b(k)t h ) 6(k) = 6(k) + 6(k)M- 1 (k)6(k), (29) 

containing the KKR matrix, Eq. (^l"[). For the reverse Fourier transform, using the right hand side of this equation, 
one obtains 

G"' =^— I d 3 k \b{k) + &(k)M- 1 (k)6(k)l e lk ' R «', (30) 

"BZ JBZ 

in which the k summation has been replaced by an integral over the Brillouin zone (BZ) with volume f^ez- 

By now, expressions for Q vold and Q have been derived for centers j of host cells only. Interestingly, the generalized 
relation 



-l 

7void. nn' 



= g™' - Q nn (t h + Q) nn Q nn ' (3i) 

turns out to be valid also, carrying arbitrary position labels n and n', which can refer to host as well as alloy positions. 
The host Green's function matrix Q nn is calculated from 

G nn ' = ^— I d 3 k[b(k,R nn ,) +6(k,R„)M- 1 (k)& T (-k,R„0]e lk - R ""' , (32) 

"BZ JBZ 

in which the matrix 6(k, R) is defined in Eq. (|l9|). The derivation is given in Appendix [b|. 

Eq. (Ill), together with Eqs. (pll) and (p5^), are the key equations of the theory presented. The alloy wave function 



coefficients Ck P L can be calculated, if the host wave function coefficients cjyfj and the matrix g vmd are known. The 
host wave function coefficients can be calculated using the KKR equation (|2^). The matrix g vold can be calculated 
after evaluation of the Brillouin zone integral in Eq. (|32|). 



III. COMPUTATIONAL DETAILS 



In this section a few quantities which are necessary in the calculation are discussed. As mentioned in Sec 
transport relaxation time r is estimated from the measur ed res istivity. The procedure is given in Subsection [II 



the 
The 

construction of the potentials is discussed in Subsection [IIIB|. Finally, the procedure for the extraction of measurable 



>cc. E| 
[II 4 ' 



quantities from the calculations is explained in Subsection III C . 



A. The transport relaxation time r 

Upon accepting the widely used isotropic transport relaxation time approximation for r, the bulk conductivity (or 
the inverse bulk resistivity) for cubic crystals is given bjo 



dS k v k . (33) 



__ 1 _ e 2 r 

p Ait^h 3 J F s 

The integral of the electron velocity v k over the Fermi surface is directly computable from the electronic structure 
of the host metal. The only quantity left is the transport relaxation time r, to which the wind force is directly 
proportional as can be seen from Eqs. (g) and (||). Combining Eqs. (|l]), (j|), (||) and (33), the temperature dependent 
wind valence can be written as 

^wind(T) = (34) 

The temperature dependence of K comes from the temperature dependence of the Fermi-Dirac function appearing in 
Eq. (^) and can be neglected up to the melting temperature of most metals. The remaining temperature dependence 
of thej-sdnd valence through the resistivity is welLJcnown from the ballistic model by Fikst3 and Huntington and 
Gronetj and from later more sophisticated modelsEZl. 

It should be noted that the calculated transport relaxation time can differ considerably from the free electron 
transport Drude relaxation time, r FE , arising from the equation 

1 ne 2 r FE , , 

a=- = . (35) 

p m 



t can be as much as a factor of 17 larger than r , as was found for palladiu 



imS 



G 



B. Potentials 



Both the host and alloy potentials are modelled by muffin-tin potentials!!^. The corresponding phase shifts, c.q. 
the t matrices t h and t p , serve as input for the computational procedure. The alloy potentials are constructed only 
for two positions of the migrating atom: the initial position and the position half-way, the saddle point. When the 
atom is somewhere else along the path, which runs along a straight line in the < 110 > direction, the phase shifts are 
calculated using the interpolation formula 

ri(s) =cos 2 M<(0) + sin 2 ( 7 r S )^(i). (36) 

The variable s runs from to 1 along the path, so rjf (0) = rjf' lmtl and halfway rjf (| ) = 7 yP> saddle P° lnt j w here p refers 
to one of the perturbed atoms in the cluster or to the migrating atom. The s dependence of the interpolation formula 
guarantees a smooth behaviour of the phase shifts at the saddle point. 

A parameter, which influences the results of the calculation, is the muffin-tin radius. The choice of this parameter is 
bound by two conflicting conditions. On the one hand it must be small enough to leave some room for the migrating 
atom, while on the other hand it must be large enough such that the potential still reproduces the electronic properties 
of the metal as much as possible. For FCC metals a muffin-tin radius of 0.325a is an optimal valueEl _ 

A real physical system is charge neutral. ThiSpGan be established using a generalized Friedel sum ruleEil and a 
shifting procedure proposed by Lasseter and Sover£3 and successfully applied in numerous systems. The potential of 
the atom is shifted by a constant value, which can be interpreted as the addition of charge on the muffin tin sphere. 
For the different configurations during the jump various choices are possible for the shift procedure. If the migrating 
atom is in its initial position, the simplest choice is a shift of the potentials of all nearest neighbours of the vacancy 
by the same value. In the saddle point configuration all nearest neighbours of the two small moon-shaped vacancies 
surrounding the migrating atom can be shifted. We want to point out that in the latter configuration the small 
vacancies are not accounted for, which means that the shifting procedure also must correct for the corresponding 
loss of charge in that region of space. This makes it unclear, what the shifting procedure means for the electronic 
properties. Due to these uncertainties and the lack of self-consistent potentials we therefore have decided to use the 
phase shifts corresponding to the unshifted potentials. Test calculations show, that the potential shift leads to a 
maximum change in the wind valence of about 10 % and in most cases of only a few per cent. In view of the limited 
accuracy of most measurements this is acceptable for the time being. In principle full-potential calculations can reveal 
detailed information of the charge state of the impurity cluster. 



C. Comparison with experiment 

The calculated wind valence is a position dependent second-rank tensor. In order to compare with experiment the 
tensor has to be reduced to a scalar. The relevant component of the force is the one in the direction of the migration 
path, which we indicate by s. Averaging over all orientations of the lattice, simulating the polycrystalline samples 
used in practice, we yield for the scalar wind valence 

Z-wind = S T ■ j?wind ■ S. (37) 

This scalar value must be averaged over the path, because the average force is the work done by it during the jump 
divided by the length of the path. By this a frequently used assumption, in which the wind force is taken as the 
average of its values at the initial and saddle point positionnJ, can be tested. 



IV. RESULTS FOR FCC METALS 



The formalism described above is applied to metals with the close-packed facc-ccntcrcd cubic (FCC) and hexagonal 
close-packed (HCP) structures, the latter being replaced by its equivalent FCC structure. Wind valences are calculated 
for diffusing atoms in aluminum, the noble metals and the 4d transition metals and are compared to experimental 
data and results from previously published computational studies. 
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A. Electromigration in aluminum 



1. Self- electromigration 

In Fig. ^| the solid line shows the variation of the scalar wind valence of a host atom in aluminum metal along 
the migration path. In addition the three partial wave contributions to Z wln <i are shown, namely the sp, pd and 
df contributions, depending on the differences r/i + i — rji as they appear in the expression for the matrix element 
< ~&k\ — VRpUpl^fc > in Eq. (||). The pd term is dominant which is consistent with the values of the phase shifts of 
the moving aluminum atom, being 0.338, 0.395, 0.051 and 0.002 for s, p, d and / respectively at the initial position 
and 0.370, 0.430, 0.056 and 0.002 at the saddle point. A sinusoidal behaviour of the wind valence along the path would 
justify the averaging procedure, which uses the values in initial and saddle point position only. Such a behaviour 
seems to be obeyed in the figure. In Table | results for the initial and saddle point position, their average value 
and the average over the path are listed. In the first three columns the three partial contributions are shown. The 
fourth contains the sum of the three. In order to obtain these values the transport relaxation time of r = 69 a.u. hsf 
been calculated using Eq. ([33]), assuming the resistivity to equal the phonon part at 800 K, namely p — 8.6 /iSlcmEJ. 
Comparing the third and fourth row one sees that the two point average Z w - md differs 10 % from the path average 
value. In all tables values for the quantity K 7 defined in Eq. (|34|), and Z w - m< i/T are given too. The advantage 
of specifying K is its weak temperature dependence (only broadening of the Fermi-Dirac distribution) , while Z w ind 
strongly depends on temperature through p. The variation of K turns out to be less than 1 % in a temperature 
range of 1000 K for aluminum, so K can be considered constant and is shown in the sixth column of Table [ij The 
theoretically purest quantity to be extracted from our calculation is Z^^/t, shown in the last column of this table. 
The calculation of this quantity requires no model for the resistivity whatsoever. A disadvantage of this quantity is 
the difficulty in comparing ii—with experimental quantities. 

Very recently Ernst et al.E3 presented experimental results for bulk electromigration in aluminum. They found an 
effective valence varying from -5 at 683 K to -3.3 at 883 K. When we use the combination of the Eqs. (Q) and ( pi[ ) 

Z*(T) = Z dilect + -?- (38) 

in order to analyse these measurements, this leads to values of K = —45 fiflcm and direct — +1-4. The value of K 
is about 1.5 times larger than the calculated value of —29 [iQcm (see Table Q), which is acceptable. It is interesting 
that the direct valence is very close to the value of half the chemical valence as predicted by Bosvieux and FriedelEi 

SorbcllccS has performed model-pseudopotential calculations for self-electromigration in Al and found a K value 
of —112 /ificm. Although a pseudopotential formulation is considered to be suitable for a nearly free-electron system 
hke-,alurninum, his value is rather different from our calculated value and the experimental value found by Ernst et 
al.Ea. On tic other hand his value agrees approximately with values given by LoddingcH, who refers to measurements 
by PenneyES. With respect to these experiments we point out, that all available theoretical models for bulk electro- 
migration lead to a temperature dependence of the effective valence according to Eq. (|38|). The measurements of 
Penney cannot be fitted by this formula. 

An interesting phenomenon is the influence of impurities on the electromigration properties of aluminum. The 
best-known example is the addition of small amounts of copper, which reduces electromigration induced damage 
effectively In the hope that other impurities have a similar but stronger effect on electromigration in aluminum, 
experimentalists tried palladium and silicon on an ad-hoc basis, but without success. In order to contribute to a 
microscopic explanation we investigated the influence oLcopper, palladium aniLsilicon impurities, that are located 
in the impurity cluster near a migrating aluminum atomo. However, we foundEl that on average the effect is small, 
although the presence of impurity atoms at particular positions can give rise to a considerable reduction of the wind 
valence. The presence of a palladium atom induces a reduction of about 10 % when averaged over all positions 
neighbouring the migration path (shaded atoms in Fig. ^). The effect of silicon and copper is even smaller. 

The presence of impurities changes the number of valence electrons in the system. Within a rigid band model 
this corresponds to a change of the Fermi energy. A vacancy or a copper atom lowers the Fermi level, while silicon 
raises €f ■ The addition of 1 %-pf copper reduces the average number of valence electrons per atom from 3.0 to 2.98. 
According to our calculations^] this reduces the wind valence only by a few percent. The conclusion is, that copper 
atoms do not reduce the wind valence of aluminum atoms directly through electronic effects. 

Finally we mention the effect of impurities due to their contribution to the electrical resistivity, which causes a 
decrease of the wind valence. An addition of 1 % of copper induces a resistivity increase of about 10 % and the wind 
valence will be reduced by about 10 %. 

We conclude that all impurity effects considered above do not induce a dramatic change in the value of the wind 
valence. Therefore the reduction of the electromigration induced damage cannot be attributed to them. 
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2. Impurities in aluminum 



Besides effects on the electromigration properties of host atoms, impurity atoms can electromigrate themselves. In 
principle, impurities can migrate in a direction opposite to host atom transport. However, the path dependent wind 
valences of copper, palladium and silicon, shown in Fig. |^, all have a negative sign and are larger in magnitude than 
■Zwind of aluminum, which is shown for comparison. This rather accelerates the vacancy transport than slowing it 
down. The average wind valences of copper, silicon and palladium are factors of 1.7, 2.8 and 7.8, respectively, larger 
than Zwind of aluminum. This is in partial contradiction with measurements, quoted by Ho and KwoktJ in their review 
article, giving copper a smaller wind valence than host aluminum. We note that the calculated ratios are independent 
of the model for the calculation of r, because t is a host quantity. The fact that the wind valences of the impurities 
have the same sign is not surprising, because the direction of the force is mainly determined by the direction of the 
motion of the charge carriers. It is worth noting, that the assumption of a sinusoidal variation of the wind valence 
along the path does not hold both for copper and palladium, which can be seen clearly in Fig. [|. This makes that in 
general the average using only the initial and saddle point position values is not a good measure for the average over 
the entire path. 



B. Electromigration in the noble metals 



In this subsection self-electromigration in the noble metals copper, silver and gold is studied as well as impurity 
migration in silver. 



1. S 'elj '-electromigration in copper, silver and gold 

The variation of K along the path for self-electromigration in the noble metals copper (solid line), silver (dotted) 
and gold (dashed) is shown in Fig. ||. A sign change from positive at the initial site to negative further down the 
path is observed. This cannot be understood in terms of electron and hole conduction and a free-electron-like model 
is not able to reproduce such behaviour. 

The calculated K values listed in Table [n] cannot be_mmpared with experiment, because as far as a temperature 
dependence of the effective valence has been measuredEirtJ, the results cannot be interpreted in terms of a constant 
K. But it is interesting to compare them to K values calculated by GuptaEJ and Sorbellcc3. Our values for copper 
and silver compare very well with the ones of Gupta, while our value for gold is a factor of two smaller. The values 
of Sorbello are very different from our results, which is not surprising in view of the pseudopotential method he used. 
However, similar trends are found in our and in Sorbello's calculations: gold shows the largest value of K, while silver 
shows the smallest. In the calculation of Gupta the K value of ppuuer is the smallest. 



In Table III we give calculated -Zwind and measured Z* valuestEfO at a given temperature T, while the last column 
contains the quantity Z w i n d/r. The calculated wind valences are of the same order p£ magnitude as the experimental 
effective valences, except for the much larger Z* of Doan and Brebec for silverEJ. In order to make a detailed 
comparison between theory and experiment more experimental information is needed. 



2. Impurities in silver 

This section is devoted to the wind valence of palladium, silver, cadmium, indium, tin and antimony impurities in 
silver. The series is part of row 5 in the. Periodic System. The calculated wind valence essentially follows the trend 
of the experimental residual resistivityEa as is shown in Fig. where both quantities are given as a function of the 
chemical valence. Both quantities are normalized with respect to the largest values, which occur for antimony. Note 
that the wind valence of the cadmium iraaurity is smaller than the one of the host atom. 

With the exception of palladium, DoarfcS investigated this series of impurities experimentally and found the effective 
valence to be linear with z(z + 1), where z — Z lmp — Z^ s , with Z- lmp and Z^ s the chemical valence of the impurity 
and silver atoms respectively. This dependence is a result of Mott's theory for Born scattering, which leads to a wind 
valence of an impurity atom proportional to z 2 at its initial position (neglecting backscattering from the vacancy) 
and to Z? mp = (z + l) 2 at the saddle point. The calculated Z w i n d and the measured Z* at about 1150 K are given as 
a function of z(z + 1) in Fig. [| Measurements and calculations show a similar trend, but the values do not agree. 
Neither the measured Z* nor the calculated Z w ind of the silver atom show the z(z + 1) dependence. 







More information about the trend in the wind valence can be extracted by considering the sp, pd and df contributions 
separately. They are plotted in Fig. ^ for the initial and saddle point position, respectively. A transport relaxation 
time of t = 243 is used, corresponding to a resistivity of p = 4.9 /ificm at 800 Kca. The pd term turns out to 
determine the overall behaviour. It develops as soon as the p states of the impurity atom are occupied, starting with 
indium. Note that for the saddle point configuration the sp term vanishes when the 5s orbit is completely filled, at 
cadmium, and leads to a relatively low wind valence of cadmium. A possible origin of it can be the following. In the 
saddle point configuration the impurity replaces two host atoms. For an impurity atom with filled 5s orbitals, two 5s 
electrons are present, just as in the configuration with host atoms in initial and final positions, which is just the host 
configuration. Such an impurity at the saddle point resembles an unperturbed host state as far as the s-electrons are 
concerned. It may be the case that scattering involving these orbitals is weaker than when just one 5s electron is 
present. 

C. Self-electromigration of the Ad transition metals 

We also applied the theory to the 4d transition metals, mainly as a challenge to experimentalists. To our knowledge 
only electrotransport in zirconium has been measuredd, showing a small positive effective valence of 0.3. A small 
positive valence has also been measured for the 5d transition metal platinum. Only metals with a close-packed 
structure are considered in this series of calculations. Metals with the hexagonal close-packed (HCP) structure, 
namely yttrium, zirconium, technetium and ruthenium, are treated as FCC metals. 

Remarkable variations of K along the path occur, as shown for rhodium and ruthenium in Fig. One observes 
minima, maxima and sign changes. Simple models will fail to describe such features. It should be stressed that the 
variations occur in spite of the smoothly varying phase shifts, so they are likely to occur due to multiple scattering 
effects. 

The K values averaged over the path are given in Table |v|. In our procedure for the construction of the potentials 
the atomic configuration is important. The electronic configuration of a free atom sometimes differs from that of an 
atom embedded in a crystal. This is the case for palladium, which has a 4d 9 5s 1 configuration in the crystal, while 
the free atom has a 4d 10 5s° configuration. The electronic configuration of a host atom in a rhodium crystal lies 
somewhere between the 4<i 8 5s 1 and 4ci 9 5s° configurations. As can be seen from the table the value of K depends only 
slightly on the configuration used. So the wind valence appears to be only weakly dependent on the precise electronic 
structure obtained from non-selfconsistent potentials. Although it remains interesting to investigate the influence of 
self-consistency we do not expect large effects for a metal like rhodium. Looking along the series in Table |y| no clear 
trend is observed in the K values. 

Also shown are the wind valences at a certain (high) temperature and the transport relaxation time r used to 
calculate it. The effective valence of zirconium, which was measured to be +0.3, is not in contradiction with the value 
of -0.4 for the wind valence However, it should be mentioned that the experiments were done at high temperatures 
in the a-phase of zirconium. The wind valence for this phase, having the BCC structure, will be calculated in the 
near future. As can be seen in the table, the wind valence in the 4c? transition metals turns out to be small. On the 
other hand the chemical valence can be rather large and therefore the effective valence is dominated by the direct 
valence and varies only slightly with temperature. Hence, these transition metals are suitable for the experimental 
determination of the direct valence, which has been the subject of a long-lasting controversy^!. 

V. CONCLUSIONS 

We have improved a Green's function method for the calculation of the electronic structure in dilute alloys. The 
formalism has been applied to the calculation of the wind force in the case of substitutional electromigration. We 
have focussed on FCC and HCP metals with the HCP metals treated in the FCC structure. 

The calculated wind-valence for self-electromigration in aluminum is in acceptable agreement with recent measure- 
ments of Ernst et alHj. We have also investigated the effect of the presence of impurities on the wind valence of 
aluminum. Neither the presence of an impurity atom near the jump path nor the impurity induced increase of the 
total resistivity of aluminum induces a dramatic change of the wind valence. So, changes in the electronic structure 
due to the presence of an impurity are not the reason for a reduction of ^ W ind- 

The wind valence for self-electromigration in the noble metals shows a sign change along the path. Such a behaviour 
cannot be reproduced by simple models. Path averaged values show qualitative agreement with experiment. 

Calculated wind valences of cadmium, indium, tin and antimony in silver roughly follow the measured residual 
resistivity and are approximately linear with z(z + 1), where z is the difference in chemical valence between the 
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impurity and the silver. Such a trend was predicted within Mott's theory for Born scattering in metals and has been 
observed experimentally, although the precise values of the wind valence do not agree very well. 

Finally the wind valences for self-electromigration in the Ad transition metals show a large variation in size and 
sign. The average values often are small due to cancellation of the wind valence along the migration path. The 
effective valence therefore is dominated by the direct valence and will depend on temperature only weakly. Hence, 
measurements on these metals are suitable to decide which model for the direct force is the correct one. 
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APPENDIX A: 



First Eq. @ will be derived from Eq. @. In Eq. (|l4|), written in the position representation, one substitutes the 
local Schrodinger equation 

^(x)^ fe (x + Rp) = (V 2 + £)<J fc (x + Rp). (Al) 
Then one applies Green's theorem and uses the equation for the void Green's function 

(V 2 + £)G void (r,r') = <5(r-r'), (A2) 
as it holds inside the void region. By that the void wave function can be written as a sum of spherical surface integrals 

n oid (* + R P )=J2j dS p' ■ (V'G void (x + R p ,x' + R p /))* fc (x' + R p ,) - 

G void (x + R p ,x' + R p ,)V'* fe (x' +Rj/). (A3) 



Using Eq. (|llj) for the void Green's function and the form of the local basis functions at the boundaries of the cell 
in Eq. (|^) Eq. ( [X^) can be elaborated for small x. This leads to the expression ([l4|) for the void wave function 
coefficients of Eq. (|10[) in terms of the alloy wave function coefficients of Eq. (|j). 

Now we want to illustrate the derivation of Eq. (|l6|), the relation between c^° ld and c^°^. This can be found by 
working out the void wave function for a position r = x + Rj = x p + R p , which is a position in the alloy cell labeled 
by p and in the host cell labeled by j, 

n° id (r) = £ c™ d jL(x p ) = £ c™ d j L (x). (A4) 

L L 

The two Bessel functions centered at different positions are related by 

Jl(x) = ji(x p + Rpj) = 3l'(Xp)Jl>l> (A5) 

L' 

in which the matrix J PJ is given by 

J P J L , - Am 1 - 1 ' ]T rCW,j£"(B«). (A6) 



Substitution of this equation in Eq. ( A4 ) automatically leads to Eq. (|l( 
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APPENDIX B: 



It will be shown how the infinite angular momentum summation in Eq. ( |l7| ) can be carried out analytically, leading 
to Eq. In addition Eq. (|25| ) is generalized to Eq. (|3l|), valid for arbitrary position labels. 

First a relation between present in the left hand side of Eq. (0), and g void >.w" will be derived. To that 

end the void Green's function G vold (r, r') for the positions r = x + Rj = x p + R p and r' = x' + Rj/ = x p , + R p /, see 
Eq. (|ll|), is rewritten using the free-space Green's function fll^). It follows that 

G void (r,r') = 

G°(r,r')+X;^(^)(^' PP ' -<0jl(x;0 = 



LU 



G°(r,r')+£ 3L^){& 3 " - ^)jl(x'). 



(Bl) 



LL' 



Applying Eq. (Bl) and using Eq. (A5) straightforwardly leads to 

^void,pp' _ qpp _|_ jpj rgvoidjj' _ gjj'\jj'p' 

Now the right hand side of Eq. ( frjf ) is rewritten using the KKR equation (^0|) and Eq. ( pi] ) 

2 J P3 (l - g void t h ) jj ' cl°f = J P3 b(k)t h c h k ° st - J p ^ voidjJ "t h c^ st . (B3) 

j' y 

The first term on the right hand side of this equation can be evaluated using the expansion 

/)"""' = J pj B»', (R pj < R jf ) (B4) 

which can be derived from an expansion for the Hankel function similar to Eq. ( |A5| ) for the Bessel function. From 
the definition of b(k) = 6(k, 0), see Eq. (|l9|), it follows that 

J pj b(k) = b(k, R p )e ik R « - B P3 . (B5) 

Substituting this equation in Eq. (B3) the first term on the right hand side of Eq. ([l8| ) appears, and one has to deal 
further with the terms 



/}'■'',),, + jwg™ d ^")t h c ^t. 



(B6) 



The sum of matrices IP-' d, r + JPig void >n' can be written as B™' + JPi(g void >n' - B"'). The latter sum is equal to 
^void, PJ ' ; given by Eq ^ (|b2|) w ith p> replaced by a host label f . By this Eq. Q o btains the form of Eq. @. 



(B7) 



Now we turn to the derivation of Eq. ( |3l| ) . Substitution of Eq. (|25|) in Eq. ( |B2| ) yields 

gvoid.pp' = B p P ' + jpj(gjf _ B"')P' V ' - J P3 G jil {t^ 1 + Q) jlh Q hj> J j ' p ' ■ 
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Regarding Eq. (BOh for Q 33 one has to _apply Eq. (B5) for a further reduction. Using the straightforward generalization 



of the host Green's function matrix (|30 



defined by 



g PP < _ 1 f dSfc^R^) + 6(k,R p )M- 1 (k)6 T (-k,R P 0]e ik ' R « 

ilBZ JBZ 

the following relations can be derived 

jpo fgjf _ B 33 ')J 3 ' P ' = Q pp ' — B pp ' — Q P3 't h B 3P ' — B P3 t h Q 3P ' 
+B P3 t h (t h ~ 1 + c,V •'/"/,"''"' 

jpjgooi = gpoi _ i ;>■'/•■ ii^ ; + gyh 

ghfjfp' = ghp' _ (t^ 1 + g) 323 ' t h B 3 ' p ' . 
Substituting these equations in Eq. (B7) gives the desired form ( |3l| ) 

gvoid.pp' _ gpp _ gpji ^ h_1 + 0) . . Q 32 P' 



(B8) 

(B9a) 
(B9b) 
(B9c) 

(BIO) 
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TABLE I. Results for self-electromigration in Al. A r value of 69 a.u., based on a resistivity p = 8.6 ^Qcm at 800 K has 
been used. 





Zwind(sp) Zwind(pd) 




Z w ind(df ) 




Z 




•^wind/T" 


initial position 
saddle point 
two point average 
path average 


0.05 -1.24 
-0.13 -6.45 
-0.04 -2.61 
0.02 -3.54 




-0.12 
0.31 
0.13 
0.17 




-1.31 -11.3 
-6.27 -54.0 
-3.79 -32.7 
-3.36 -28.9 


-0.019 
-0.091 
-0.055 
-0.051 


TABLE II. 


Calculated K values for the noble metals compared with res 


ul 


ts by GuptaEl and SorbellcR 




present work 






Gupta 


i>:i 


SorbellcEf 


Cu 
Ag 
Au 


-33.3 
-24.7 
-42.6 






-31.5 
-33.1 
-83.2 




-196 
-171 
-229 


TABLE III. Wind and effective valences for the noble metals. The 2 nd and 3 rd columns contain the calculated wind valence 
and the measured effective valence, respectively, at the temperature shown in the 4 th column. The last column contains the 
quantity Z w i n( j (total) /r in atomic units. 




■^wind 










T(K) 


Zwind/ T 


Cu 
Ag 

Au 


-3.5 
-3.3 

-3.4 


-4.3SJ 
-19.sk. 

-5.l£ 

-6.eE 


J 






1300 
1150 

1289 


-0.028 
-0.021 

-0.037 



TABLE IV. The calculated K and Z w i n d at a temperature T for the 4d transition metals. The transport relaxation time r 
at that T is given in atomic units. 



K(fincm) T(K) Z wind r(a.u.) 

Y 93^2 1700 OA 25.4 

Zr -52.1 1700 -0.4 13.7 

Tc 34.8 1700 0.5 12.4 

Ru 7.0 1500 0.2 10.6 

Rh^d^s 1 ) -36.9 1700 -1.0 13.0 

Rh(4d 9 5s°) -38.3 1700 -1.0 12.7 

Pd -58.6 1700 -1.3 50.4 
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FIG. 1. The systems used in the description of the dilute alloy (left): the void system (center) and the host system (right). 
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FIG. 2. The variation of the wind valence of a migrating host atom in aluminum at 800 K. The sp, pd and df contributions 
are given separately. 
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FIG. 3. The position dependent wind valence of copper, palladium and silicon impurities in Al. 
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FIG. 6. The solid circles give the calculated wind valence of the atoms Ag, Cd, In, Sn and Sb in Ag as a function of z(z + 1), 
where z — Z at0 m — ^A g - The open circles are effective valences, measured by DoanEJ at temperatures of approximately 1150 K. 
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FIG. 7. sp, pd and df contributions to the wind valence of 4d and 5sp impurities in Ag at T — 800K. 
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